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Abstract 

We elucidate the effects of chiral quenched disorder on the scaling properties of pure systems by 
considering a reduced model that is a variant of the quenched disordered cubic anisotropic 0{N) 
model near its second order phase transition. A generic short-ranged Gaussian disorder distribution 
is considered. For distributions not invariant under spatial inversion ( hence chiral), the scaling 
exponents are found to depend continuously on a model parameter that describes the extent of 
inversion symmetry breaking. Experimental and phenomenological implications of our results are 
discussed. 
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I. INTRODUCTION 



The large-scale, macroscopic effects of disorder in statistical mechanics models and con- 
densed matter systems have been a subject of intense study for a long time now. Very 
generally, depending upon time-scales, there can be two kinds of structural disorders that 
can exist in a system, namely, annealed and quenched disorders. In a system with annealed 
disorder the impurities can diffuse freely in a system until they come to a thermal equilib- 
rium state. The time-scale of the annealed disorder dynamics is comparable or shorter than 
that for the physical degrees of freedom of the corresponding pure system, and hence the 
dynamics of the disorder becomes important. In contrast, for a system with a quenched 
disorder the impurities are fixed in particular configurations and do not evolve in time, 
and, as a result, the disorder configuration is not in thermodynamic equilibrium. Studies 
on the effects of random quenched disorder on pure systems, e.g., the classical 0{N) spin 
model 1-3| and self-avoiding walks on random lattices j4-8| are important, particularly be- 
cause of the modifications in the critical behaviour brought about by the presence of such 
impurities. Whether or not quenched disorder changes the universal scaling properties of 
the pure system is given by the heuristic arguments of the well-known Harris criterion j^. 
^r^aussian-distributed random impurities, both spatially short-ranged joj] or long-ranged 



13 



variances have been considered. For instance, perturbative renormalization group 
(RG) calculations on 0{N) symmetric models yield that for short-ranged disorder, the scal- 
ing exponents, even when they are affected by the disorder, do not generally depend upon 
the strength of the disorder 9|. In addition, experiments on liquid crystal systems reveal 
important (i.e., relevant in a RG sense) effects of quenched disorder [lO|. These studies typ- 
ically considered disordered media with variances invariant under spatial parity inversion 
(achiral). In contrast, the influence of chiral disordered media on the scaling properties of 
pure systems are yet to be explored systematically. Chirality is known to affect universal 



properties of systems; see, e.g., Refs.jll 



121] where effects of chirality is discussed on a class 



of driven nonequilibrium models. Furthermore, it is now possible to fabricate chiral liquid 
crystals by embedding nematic liquid crystals into porous chiral film made of deposition 
of helical columns of MgF2 on glass substrates. This is believed to be of use in practical 



appHcatio,., e.g., optical =w.tcM„g M- view of the=e examples and consMering the fact 

that chiral disordered substance are now prepared experimentally (e.g., chiral aerogels |14j|). 
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it is pertinent and relevant to theoretically investigate the effects of chiral disordered media 
on the scaling properties of pure systems in a general setting. 

In this article, we propose a reduced minimal model with quenched chiral disorder and 
analyse it to study the generic effects of such a kind of disorder on the critical scaling of 
pure systems. While our intention is not to model any specific real system in details, our 
results should have broad implication in understanding of media with orientational disorder, 
e.g., the effects of chiral porous media [l^ on the nematic-to-smectic A (N-A) transition 
of liquid crystals or on the scaling properties of smectic-A liquid crystals. Effects of usual 
(achiral) short-ranged disorder on N-A transitions and smectic A are both theoretically and 
experimentally well-studied 15| and a variety of results are obtained including possible de- 
struction of translational (smectic) order, but stabilisation of the smectic Bragg glass phase. 
There are other systems where effects of chiral quenched disorder is likely to be of signifi- 
cance, e.g., disordered cholesteric hquid crystals 15|, smectic A to C (A-C) transition [igI 
and superfiuids in aerogels The first one would be particularly intriguing, due to the 
fact that pure cholesteric liquid crystals are themselves chiral. However, the ensuing al- 
gebraic manipulations for the specific systems mentioned above in conjunction with chiral 
quenched disorder and its coupling with pure system variables, will be quite challenging. 
Thus, studies on a simpler reduced minimal model should be welcome. Such approaches are 
useful provided it allows one to address questions of basic principles, which in the present 
case is the effects of chirality or handedness of the disorder medium. To this effect, in this 
article we consider a variant of the usual quenched disorder 0{N) model (A^ > 2) with cubic 
anisotropy in the presence of a disorder distribution that breaks the invariance under spatial 
reflection or parity inversion and study its scaling properties near the critical point. Our 
model is essentially a generalisation of the random cubic anisotropic model introduced in 
see also Ref. |l9j for related discussions on a problem of unconventional super- 



Ref. 



conductors with quenched impurities. The lack of parity inversion symmetry in our model 
represents chirality or handedness of the disordered medium. Possible macroscopic effects of 
the chirality of the disordered medium, e.g., effects on universal scaling properties have not 
been considered before; nevertheless, formally their existence cannot be ruled out on any 
general (microscopic) ground, for they represent lack of reflection invariance of the impurity 
distribution at small scales. By systematically using a perturbative RG framework together 



with the replica formalism 



20 



21| . we find (a) for order parameter component > 4, 
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the disorder is irrelevant (in a RG sense, see below) with the second order phase transition 
in the model being described by the corresponding pure system scaling exponents, (b) for 
< 4 when the disorder distribution is reflection invariant, the system is described by a 
random isotropic (Rl) flxed point (FP), at which disorder is relevant (in a RG sense) and 
which is identical to the scalingbehaviour of the isotropic 0{N) symmetric model with 
short-ranged quenched disorder j9|, in agreement with Ref. Ig], and (c) one may vary this 
random isotropic FP continuously by tuning a non-negative model parameter A^x (see be- 
low), that describes the strength of the parity-breaking parts in the variance of the disorder 
distribution from zero. The scaling exponents that characterise the second-order phase tran- 
sition at this FP vary continuously with A^x- At the technical level, as we shall see below, 
A^x appears as a marginal operator in the ensuing disorder-averaged state of the system. In 
addition, for a non-zero A^x, near the critical point the model displays a diverging correla- 
tion length ^ and a fluctuation-corrected renormalized that depend on A^x , suggesting the 
possibility of A^x -dependent scaling near renormalized critical point in different realizations 
of the model system. Our results should be directly testable in disordered classical systems 



with carefully chosen distributions for the disorder at dimension d = 3 . The remaining 
part of the article is organised as follows: In SecHIl we propose and deflne our model. Then 
we extract the critical scaling exponents from the model in Sec JIIII In Sec JIVt we summarise 
and discuss our results. 



II. THE MODEL 

Let us begin by discussing possible requirements of a simple model that will be able to 
capture the effects of chiral quenched disordered media. First of all, if the chirality of the 
disordered medium has any discernible effect on the system, we expect two-point correlation 
functions of the physical variables that describe the corresponding pure system to display 
chirality (i.e., have parts that are odd parity under spatial inversion). This evidently requires 
two or more variables describing the disorder in the system, since only a cross-correlation 
function of two different variables can have a part that is odd under spatial inversion. 
Secondly, for an impurity distribution that displays lack of symmetry under inversion, it 
is expected to be described by more than one (frozen) fleld. Keeping these features in 
mind, and considering systems that are purely achiral in the absence of any disorder, we 
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consider a disordered version of the well-known pure cubic anisotropic 0{N) model [23j as a 
reduced minimal model for chiral quenched disorder systems. The corresponding free energy 
functional is J-" = J d'^xf with 

^ 1 

i=l 

where tq = (T — Tc)/Tc with T and referring to the temperature and the mean-field 
critical temperature, respectively. Stochastic function rj(x) represents the coupling with 
the disorder, such that = + rj(x) is the local fluctuating critical temperature for 
(pi. Further, m > and f > are the bare nonlinear coupling constants in the model. In 
the limit of v = and for all rj(x) = f(x), the microscopic rotational invariance in the 
order parameter space is restored and we get back the usual 0{N) model with quenched 
disorder. On the other hand, if all = and f 7^ 0, it reduces to the well-known cubic 
anisotropic model Furthermore, with all rj(x) = f(x) and both u, f 7^ 0, it is identical 
to the model of Ref We have kept the gradient term of the free energy functional 

[1] spatially isotropic for simplicity. This, though admittedly an idealisation, simplifies the 



ensuing algebra considerably without destroying the effects of chimlity of the disorder 24 1. 
The cubic anisotropic terms reflect any possible breakdown of rotational invariance of the 
order parameter. The relative relevance (in a RG sense) between u and v has been addressed 



by perturbative RG calculations, see, e.g., |23|: For < 4, the scaling properties are 



described by the stable isotropic FP f = 0, m/j > 0, with suffix R referring to renormalized 
quantities, and the associated scaling exponents at the critical points are identical to that 



for the usual 0{N) model (with bare v = 0) [23[- In contrast, for > 4, the cubic 
anisotropy is a relevant perturbation on the isotropic FP and the system cross over to the 
cubic anisotropic FP with Ufi > 0, vr > 0. In addition, the system can display fluctuation 
induced first- order transition due to the Coleman- Weinberg mechanism for a range of (bare) 



values of the coupling constants 25||. With this background and in order to include the 



effects of spatial reflection symmetry breaking in the disorder distribution, we assume the 
fluctuations in T^, rj(x) in our model to be Gaussian distributed with variances 

(r,(x)r,(0)),, = 2Z}5,/(x) + a,,[D5(x) + D(x)], (2) 

with the matrix aij having a structure of the form an = = 022 and ai2 = 1 = 021; function 
D{pc) = — D(— x), D is a positive constant and £) is a constant that can be both positive 



and negative. The symbol {■■■)av represents averages over the chosen Gaussian disorder 
distribution. A non-zero -D(x) thus introduces breakdown of parity in the system. Since 
D(x) is an odd function of x, its Fourier transform must be imaginary and odd in Fourier 
wavevector q, such that the Fourier transform of -D(x) is given by iD^{q), where -Dx((l) is 
a real odd function of q: -Dx(q) = ~-Dx(— q). In order for -D(x) to be equally relevant (in 
a RG sense) with D6{'x) and Z)(5(x), i.e., for them to have the same physical dimensions, 
-Dx(q) should have the same g-dependence (in a power counting sense) as D or D. Since, 
D and D are just constants, independent of q, -Dx(q) must not depend on the magnitude of 
q. Therefore, we further define an amplitude = Dy^{q)D^{q) and a dimensionless ratio 
A^x = (-Dx/-D)^ > 0. Clearly, -D(x) is an odd function of x that has the same dimension 
l/L*^ (L is a length) as the 5-function 5°'(x) in d-dimension. It is easy to work out an explicit 
representation of -D(x) in Id: Writing D{x) as the inverse Fourier Transform of Dx{q) where 
g is a one- dimensional Fourier wavevector, we find 

POO 

D{x) = i I dkDy; [exp{ikx) — exp{—ikx)] ~ D^/x. (3) 
Jo 

Although 1/x has a range longer than 6{x), it has the same physical dimension and hence 
scales the same way as 6{x) under rescaling of space, i.e., x. This paves the way for different 
elements in ([2]) to compete in an RG sense. The d-dimensional analogue of the Id form of 
D{x) above is rather complicated. Nevertheless, it should generally be of the form on 
general dimensional ground in one hemisphere, with a change in sign in the other hemisphere. 
We then ask: What are the scaling properties of our model? Our heuristic arguments and 
detailed calculations below reveal that A^x may be varied to tune the emerging critical scaling 
behaviour continuously. Before we discuss our results in details below, a few words about the 
interpretation of the structure of our model in the context of possible physically realizable 
examples are in order: The mixing of order parameter indices and spatial dependence in 
([2]), although not allowed in the usual spin models, is consistent with the Frank free energy 
for nematic liquid crystals [24] where the director field is defined in the coordinate space. 
A chiral disordered material has both positional and orientational disorder. For instance, 
in a chiral aerogel the positional disorder is related to the local density fluctuations of the 
aerogel pores, where as the orientational disorder reflects the randomness in the orientation 
of the pores. The latter one, represented formally by a quenched vector field 15|, may in 
general have a parity breaking variance, and should couple with the fields, e.g., the nematic 
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director field in the N-A or the displacement fields in the A-C transitions. While detailed 
form of such couplings are model-specific and are not necessarily as simple as we have in 
([2]), these, in-principle, should generate disorder distributions that breaks symmetry under 
inversion. A simple choice as ([2]), despite its limitations, suffices for our purposes here. In 
general, the distributions of the two types of disorders may not have any simple relation as 
they may occur independently. 



III. CRITICAL SCALING EXPONENTS 

Our model ([2]) without any disorder (all = identically) displays standard order- 
disorder transition through a second order critical point (in addition to a fiuctuation induced 
first order transition) as discussed above. To get an idea about the possible macroscopic 
effects of disorder on the pure system properties, it is instructive to first consider the predic- 
tion of the Harris criteria jsj for the present model. In order to retain the effects of the parity 
breaking part of the disorder variance while constructing the Harris criterion, we formulate 
it as given below. While here we closely follow the derivation of the standard Harris criterion 

n 

and the notations as in Ref. [3i], we nevertheless rephrase the details here again for the sake 
of completeness. To this end, we divide the system into subsystems of linear dimension ^, 
where ^ is the correlation length at that temperature of the corresponding pure system. The 
idea is to find out if the variation of the critical temperature of these regions of size ~ ^ 
becomes negligible as T — )■ Tc. Since the spins are expected to be correlated and on average 
aligned for up to a distance ~ C,, the transition temperature for the ith component (pi of 
the order parameter field of a region of size ^ may be defined as the average of T^^(x) over 
that region. We define reduced temperature t = (T — T^/T^ and local reduced temperature 
tj(x) = fi{'x)/Tc, fj(x) = tq + rj(x) for the i-th component of the order parameter field, we 
have {ti{'x))av = t. Further, as defined in Ref. jsl 

Uv = ^, I d'^xtii^) (4) 



V . 

is the effective reduced temperature of a region V = Note that we have formally allowed 
an effective reduced temperature tiv for the i-th component (pi of the order parameter 
field. This is consistent with the fact that our model allows for order parameter component 
dependent effective reduced temperature in the free energy functional [1] The variance Ajj 
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of Uv is defined as 



I av 



^ij — {tivtjv) 



c 

av 



= i^i^yJ^d'^xj|^A^?.,(x-y), (5) 

where (7jj(x — y) = {ri{x.)rj{y))av, as given in Eq. ([2]). Here, a superscript c refers to the 
connected part of the variance. Thus 



A.,- = j^d'x j^d'y(2D5i,5{^-y) + a,,[b5{^-y) + b{^-y)] 



(6) 



It is clear from ([6]) that the contribution from the odd parity part of the disorder variance 
to the variance Ajj vanishes owing to the odd parity of -D(x). In order to capture the effect 
of the parity breaking part of the disorder variance (i.e., non-zero -Dx), we modify Eq. ([6]) 
to 

A,, = 1^1^ J^d'x J^d''y(2D5,,Si^-y) + a,,[bSi^-y) + \Di^-y)\]) . (7) 

Alternatively, one may restrict the domain of integrations above to hemisphere having a 
single signature of the parity breaking part in the disorder variance. While the above 
modification, in terms of considering the absolute value of i5(x — y), is admittedly apriori 
designed to capture non-zero contributions from the parity breaking part of the disorder 
distribution, this does not alter the power counting in the integral in (|6]). For = this 
modification still leads to the well-known Harris criterion js, 5|. With our modification, 
therefore, Ajj is expected to have a contribution proportional to \D^ \ or D"^. From ([7j) we 
note that, Aij will have a part ~ ^"'^ coming from D and -D in ([2]). In addition, there should 
a part ~ .^^'^In^ (this may be shown explicitly in Id with D{x) ~ l/x) coming from the 
parity breaking part in ([2]). Proceeding as in Refs. 0, we then conclude that disorder 
is relevant as long as the specific heat exponent of the corresponding pure system a > 0. 
This condition is same as the usual Harris criteria Consider now the fact that the 

borderline of relative relevance (in a RG sense) between the short-ranged and long-ranged 
disorder (with a variance Ixl*^) is determined by the condition a = d ([3]), which yields a 
logarithmic contribution to the analogue of in Ref. js]. Thus looking at the logarithmic 
dependence associated with the variance Ajj above, it appears that the present model is at 
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the borderline between ((5-correlated) short-ranged and long-ranged disorder. Furthermore, 
for relevant long-ranged disorder the scaling exponents depend explicitly on a. Since we can 



write Ixl" 



d+S 



5\x\'^\nx with a = d — 6, 5 — )■ 0, 5 should appear as a control parameter 
in the scaling exponents. Hence, drawing on the analogy between Ajj as above and the 
corresponding expression in Ref. j^, and comparing with ([2]), amplitude (equivalently 
A^x) should appear as a tuning parameter. Thus, any correction to the critical exponents 
due to Dx must be at least 0(Dx)^ [and hence O(A^x)], since our perturbative calculations 
given below should be analytic in Dx (or A^x), where appears as an expansion parameter. 
While our arguments above are of heuristic nature and do not constitute a rigorous proof, 
they are indicative of non-trivial behaviour with finite D^] our detailed RG calculations 
below confirm this qualitative physical picture. 

In order to systematically investigate the properties of systems with quenched disorder it 
is required to average the free energy over the disorder distribution. This can be conveniently 
done using the well-known replica method 2l|] We start with the partition function for the 
free energy functional ([1]) 

N 

Z = l[Tr^^exp[-pj^{^,}]. (8) 

1=1 

Then the free energy averaged over the disorder distribution can be written as 



F = -{\nZ), 



lim 

m— i>0 



lim 

m— >-0 
fN 



m 



m 



(9) 



Here, a = 1, 2, m are the replica indices and {0"} represents m replications of the order 
parameters 0j. The corresponding m- replicated disorder averaged partition function is given 
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by (we set ksT = 1, where /c^ is the Boltzmann constant) 

m TV „ 

{z^av = n n / ^'^^ 



a=l i=l 



exp 



X exp 



X exp 



N m 



(Fx 



/, m N N m 

N m p N m 




j=l a,/3=l 



/.2 



i,j=l a, 13=1 



X exp 



N 



J2 E 0L(a:)^(a;-a:')0?/3(^') 



i,j=l a, 13=1 
i<j a^P 



(10) 



Here a, /3 = 1, 2, n are the rephca indices. Nonhnear terms with couphng constants D 
and D involve fields with different replica indices; thus, these terms lead to mixing of replica 
indices. These terms in Eq. (fTOj) originate due to the averaging over the disorder distribution. 
Our purpose is to calculate the scaling exponents r] (anomalous dimension) and v (correlation 
length exponent) near the second order phase transition, which are formally defined through 
the relation (0j(r)0j(O)) ~ ffji^ / C) where ^ ~ |T — Tc\~" is the correlation length; 

ffj is a dimensionless scaling function and averages (..) are thermal averages to be obtained 
from the disorder averaged free energy functional. Anomalous dimension r] describes the 
spatial scaling of the correlation function at T^, where as v describes the divergence of ^ as T 
approaches renormalized T^. If we take only the terms quadratic in 0" in Eq. fllOp . the model 
can be solved exactly. When the nonlinear terms are present, exact solutions are ruled out. 
The RG framework provides a systematic method to extract scaling behaviour in a nonlinear 



theory^ 



Ref. 



25 



he detailed methods of the RG framework is well-documented in literature, see, e.g.. 



261] ; see also Ref. 9| for applications of RG in a disordered system. The presence of 



the nonlinear terms requires expanding the different vertex or correlation functions in powers 
of the coupling constants. Near the critical point the perturbative corrections diverge leading 
to failure of naive perturbation expansion. The perturbative corrections are represented 



by using the standard Feynmann diagrams 



We use a minimal subtraction scheme 



together with the dimensional regularisation scheme to evaluate the diagrams, where only the 
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diverging parts of the diagrammatic corrections are obtained in inverse powers of e = dc — d. 
Here dc is the upper critical dimension, at which the relevant coupling constants become 



dimensionless. For the present model, d^ = 4: for all the nonlinearities [27|. Thus all of them 
are equally relevant. We define the renormalized coupling constants via the renormalization 
Z-factors ur = Zy^ufi~''Aa, vn = Z^vfi^'^A^, -D_r, = Z^D^^'^Ad, Dr = Zj^D^^'^Ad, where 
/i is an arbitrary momentum scale and A^ = ■^1^1/2- The Z-factors, defined to absorb the 
divergences, are used to obtain the RG equation for the correlation or vertex functions, 
which yield the scaling exponents at the different RG FPs. The FPs are formally given 
by the zeros of the RG /3-functions (3a = l^-g^ci^ ci = u,v, D, D in the present model. The 
/9-functions are found to be (there are no corrections to up to two loop orders) 

/3u = [-e + A8uR + 8ur{N + 2) + 48vr - 32Dr 

- 32Dr]ur, 

13^ = [-e + 96uR + 72vR ^ 

vr 

+ ^ - 9QDr]vr, 

/3d = [-e + 48{uR + VR)-64DR 
+ 8iN-l)^]DR, 

^R 

(3f) = [-e + A8iuR + VR)-32DR + 32^^^ 

Dr 

+ 16(iV - 2)ur - IQDr - IQ^]Dr. (11) 

Dr 

The RG FPs are given by the zeros of the /3-functions ffTTl) . It is useful to first consider the 
case with Dy^ = 0. The FPs are 

• Gaussian FP: ur = 0, vr = 0, Dr = 0, Dr = 0, 

• Heisenberg FP: ur = g^^, vr = 0,Dr = 0, Dr = 0, 

• Cubic anisotropic (CA) FP: ur = ^,vr = ''^^'^^ , Dr = 0,Dr = 0, 

• Random isotropic (RI) FP: ur = ^2{n-i) ^ = 0,Dr = Dr = 54 (]v-\') ' 

• Random cubic (RG) FP: ur = ^^f^z^^, vr = ^j^E^, Dr = ^j^, Dr = 2Dr = 

e A-N 
96 N-2- 
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Among the coupling constants, only Dji must be non-negative in order to be physically 



meaningful. Thus, both RI and RC FPs may exist only for 4 > [28|]. For > 4 the 
disorder will not be relevant any more and the system will be described by the pure system 
FPs. Although even with Dy^ = microscopically our model here is different from (and 
a slight generalisation of) that in Ref. 18 1 in having two parameters D and D denoting 
variances of the disorder distributions, note that at the RI FP, = 2Dii. This holds at 
the RC FP as well, although the overall system is no longer isotropic (due to a non-zero vr) 
in the order parameter space. Thus, generically, Dr = 20^ at the FPs that depend upon the 
disorder. This shows the redundance of having two independent parameters D and £) in a 
RG sense, and demonstrates the robustness of the results of Ref. [l8|. We, however, shall see 
below that independent D and Dr are required when Dy. ^ 0. Unsurprisingly, our results 
on the FP values of the coupling constants at the RC FP above (with -D(x) = 0) match 



unctions (fTT]) exactly 



with those of Ref. 0. In fact, with Dr = 2Dr and L'(x) = 0, the 13- 
correspond to the recursion relations for the coupling constants in Ref. [18| . At the Gaussian, 
Heisenberg and CA FPs, the exponents are well-known |^. With the RG FPs available, the 
different scaling exponents may now be obtained by using standard procedures: One begins 
by calculating the the RG Z-factors Z^. and Z^, defined via renormalised tr = Z^vq and 
(f)R = Zfp(f). The Z-factors Zj. and Z^ are obtained from the one- and two-loop Feynmann 
diagrams by using a minimal subtraction scheme together with dimensional regularisation. 
The corresponding Wilson flow functions 7,. = In Zr and 7,^ = In then immediately yield 



the exponents u and t] respectively. For further details we refer the reader to Refs. 25|, |26 |. 



The scaling exponents at the RI FP are modified by the quenched disorder and are given by 
_ {N + 2)e^ {A-N){N + 2)e^ {N - Afe^ 

^ ~ 00/ AT TT9 TJTTlvr TT9 



32(Ar-l)2 64(A^-1)2 256(Ar-l)2' 

1^2- , (12) 

8(A^-l) ^ ^ 

These are identical to those at the random FP in Ref. [3]. The remaining critical exponents 
may be obtained from the exponent expressions ([T2|) . We perform linear stability analyses 
around the RI and RC FPs; we consider only A^ < 4, for which these FPs are physically 
meaningful [29 1. Since 2Dr = Dr is maintained at both RC and RI FPs, we are concerned 
with linear stability in the {ur, vr, Dr) space. Further, we concern ourselves with second 
order phase transitions only. First, the RI FP: We find for the eigenvalues of the linear 
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stability matrix 



A 



(4 - A^)e (4 - A^)e 



(13) 



' 4(A^- 1)' 4(A^- 1)' 

Thus, for A^ < 4 all the eigenvalues are positive, and hence stable in all the three directions. 
The discussions of the stability of the RC FP in Ref. 18|] directly apply here as well. For 
the sake of completeness, we show the eigenvalues up to 0(e) at the RC FP 

(A^-4)e (4-A^)e 



A 



(14) 



6(A^-2)' 6(A^-2)' 

Thus, two of them diverge at A^ = 2. It is thus gene rally expected that flow lines from near 
the RC FP flow to the (stable) RI FP. See Ref. 18|] for detailed discussions on this. 

Having known the FP values of the coupling constants for D-^ = or A^x = 0, we 
now obtain the corresponding values when Ny_ > 0. For simplicity of the ensuing algebraic 
manipulations, we find corrections to the (A^x = 0) FP values of the coupling constants up to 
0{Nx), assuming a small A^x • Thus the finite A^x FP values, as we write down below, are not 
going to be quantitatively accurate for A^x ~ 0(1). For = 0, the zeros of the /3-functions 
(3d and as given in Eqs. (ITT!) , yield Dr = 2Dr generically at the RG FPs. This then, 
in turn, leads to Dr = = Dr as a FP solution, regardless of its linear stability properties. 
Thus, pure system properties may be restored in the large scale provided Dr = = Dr 
is a stable FP. However, when D^ 7^ 0, two observations may be made immediately: (i) 
Dr 7^ 2Dr and (ii) Dr = and Dr = are no longer solutions of = and f3£) = 0. 
Thus, pure system behaviour is not expected to be observed. We obtain modifications to 
both RI and RC FPs up to O(A^x) separately. First the modified RI FP: After some straight 
forward algebra, we find 



u 



R 



'R 



32(A^- 1) 



2(4 + AT) 



:(4 - N)N^ 



Dr 
Dr 



16(4 + A^)(A^- 1 
e(4 - A^) 



128(A^- 1) 
e(4 - A^) 
64(A^ - 1) 



1 + 



2(4 + A^)' 

(A^ + m-x 

2(4 + A^) 



(15) 
(16) 
(17) 
(18) 



Here, A^x 



El 



> is a dimensionless parameter. Notice that: (i) vr picks up a small 
non-zero value O(A^x) and (ii) Dr 7^ 2Dr. Thus the modified RI is no longer isotropic; 
however, the departure from isotropicity is small, O(A^x)- It is clear from Eqs. (llStilSp that 
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for 1 < < 4 the FP values are physically meaningful, e.g., Dr is positive definite, for any 
values of (small) A^x; we, therefore, consider only the range 1 < < 4 below. Corrections 
to 0{Ny) for the RC FP may also be obtained, which we do not show here explicitly; see 
Ref. jisj for the details in this context. For small A^x and 2 < A^ < 4 this modified RC 
FP remains unstable. The general picture of the second order phase transition remains 
unchanged at such finite but small A^x- When A^x ~ 0(1), the FPs are expected to be 
substantially modified, which we do not discuss here. We are unable to comment about the 
physical picture for A^ > 4 on the basis of the FPs given by fllStllSp , since Dr < for A^ > 4 
from Eq. ([T71). 

The critical exponents rj and z/, which may be evaluated following the procedure outlined 
above, at the modified RI FP are: 

_ {N + 2)e'' (4- Ar)(Ar + 2)e' [N-Afe'' 

^ ~ or./ AT TvJ FTTaT TT^ 



32(A^- 1)2 64(A^-1)2 256(Ar-l) 

(4 - N^e^ 



256(A^- i; 



ATx, (19) 



1 3A^e 3(4-A^)e,, , , 

As usual, the critical exponent expressions (fT9l) and ( 120|) yield all other critical exponents. 
The scaling exponents given in ( IT9l) and (1201) constitute the principal results of this work. 
From Eq. (fT9|) . for A^x = 0, is always positive (we consider 1 < A^ < 4). For small A^x 7^ 0, 
rj remains positive, but its value decreases. Thus the correlation function for (pi decays more 
rapidly with spatial separation for A^x = than for A^x 7^ 0. It apparently suggests that rj 
can be brought to zero for sufficiently large A^x; however, since our expression ( IT9l) are valid 
for small A^x only, we cannot say anything conclusively when A^x is large ~ 0(1). Similarly, 
from Eq. ( l20i) . z/~^ remains positive for small value of A^x and it has a larger value when 
A^x > than for A^x = 0. Since u is an explicit function of N^, the measured diverging 
correlation length ^ ~ \T — Td^" will depend on the value of Ny.; ^ for A^x > is larger than 
that for A^x = 0. Furthermore, notice that u picks up an A^x-dependent part at 0(e), where 
as the same for t] appears at O(e^). Thus, experimental detection of any A^x dependence will 
be revealed much more clearly in measurements of the correlation length near renormalized 
T 

What could be an upper bound of A^x? By demanding positivity of the eigenvalues 
of the disorder variance matrix, we can enforce a bound on (or, equivalently, an upper 
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limit on A^x)- Eigenvalues Ag in general depend both on A^, the number of order parameter 
components, and D, D^- We take = 2 as a specific example. The variance matrix M 
of the disorder distribution is given by 

2D D + iD^ , 
M= . M • (21) 

D - iDy 2D 



The corresponding eigenvalues are Ag = 2D ± y (Z)^ + D'^) > 0. This gives a bound on the 
off-diagonal elements of the variance matrix M. This is consistent with the limit Dx = 0, 
where Dr = 2Dr at the FP. Similar exercises may be under taken for higher values of A^, 
which we do not discuss here. Finally, our claim of continuously varying universal properties 
rests on the possible continuous variation of Dy., and hence of A^x • We have shown, in our low 
order (two-loop) renormalized perturbation theory, that there are no fluctuation corrections 
to Dx- We believe this holds to any order in the perturbative expansion. To see this 
notice that any non-zero fluctuation corrections to -Dx(q) niust be an odd function of the 
its wavevector argument. In order to have such a non-zero correction one must have an 
odd number of -Dx(ct) vertex in the diagram. Since all internal wavevectors are integrated 
over, such a contribution will vanish in the limit of vanishing external wavevector. Thus 
any putative diagrammatic corrections to -Dx(q) vanishes and hence A^x should appear as a 
dimensionless marginal operator to any order in perturbation. This is technically similar to 



a marginal operator that exists in the models of Refs. [ll|, |12| . Numerical verification of our 
results on equivalent lattice-gas models requires generation of A^ stochastic functions having 
variances as given by ([2]). This may be conveniently done by following the method outlined 



in Ref . 



12| (see also Ref. [30| for a general discussion on related issues) 



IV. SUMMARY AND OUTLOOK 

In summary, thus, we have proposed and studied a variant of the classical cubic anisotrpic 
0{N) model with short ranged quenched disorder having a parity breaking part, with a 
strength parametrised by D^ (equivalently by A^x), in its variance, as a reduced minimal 
model to study the effects of quenched chiral disorder on the scaling properties of pure 



systems. For our work, we use a generalisation of the model used in Ref. 



m 



The truly 



novel result from our reduced model is that the explicit and continuous dependence of 
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the scaling exponents on A^x 7^ 0, with our results reducing to those of Ref. jl8| in the 
limit Nx = 0. These A^x-dependent results on the scaling exponents are qualitatively new 
contributions made in the present work. We have worked out the dependences of the scaling 
exponents up to O(A^x) for simplicity using two- loop e-expansions. If the disorder is long- 
ranged j^, i.e., if the magnitudes of D, D, Dy;(k) have k-dependent parts ~ k^^, y > 0, 
in addition to constant parts, preliminary calculations analogous to that here reveal that 
a parameter analogue to A^x will appear in the scaling exponents. Thus, the effects of 
a parity breaking part appears to be quite robust. Assuming the disorder in our model 
results from a variety of microscopic sources, the tuning parameter A^x may be interpreted 
as a measure of the relative concentration of the microscopic impurities that causes a non- 
zero parity breaking variance in the resulting impurity distribution. Hence, our results are 
illustrations of dependences of the scaling exponents on the relative concentration. Our 
reduced model and the results that follow are not directly applicable to any specific system 
we discussed (e.g., the N-A transition), due to the underlying simplicity and idealisation of 
the model. Nevertheless, the broad qualitative picture that chirality of impurity distributions 
may be relevant in determining the universal properties of disordered systems is sufficiently 
general and expected to be observed in more realistic and specific experimentally accessible 
systems. From a general point of view, our results open up the possibility of a new paradigm 
in the scaling properties of quenched chiral disorder systems. Notice that our results do 
not contradict the recent results that critical exponents in the disordered Ising model are 
independent of impurity concentration along the transition line between paramagnetic and 
ferromagnetic phases 311], since = (hence A^x = 0) identically for the Ising model 
with A^ = 1. Many experiments tend to suggest smearing of phase transition in a variety 
of disorder systems j32:]. Our results on A^x-dependent correlation length as T —> Tc and, 
similarly, A^x -dependent renormalized are reminiscent of a system with a continuous 
spectrum of critical points and relevant diverging length scales (near Tc) and critical points, 
and thus loosely resemble a smeared transition (although there are no formal connections). 
Measurements on model systems (for A^ > 2) with different realizations of the disorder 
having different A^x? corresponding to a given impurity distribution, naturally leading to 
a broadening of the measured values of the scaling exponents. Although at present no 
experiments on chiral disordered systems are available to our knowledge, we expect such 
experiments (e.g., N-A or A-C transitions) should be performed in the near future by using, 
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say, chiral aerogels [ij]. Numerical simulations of our reduced model with appropriately 
chosen disorder distributions should be useful. From the point of view of the notion of 
universality in both equilibrium and nonequilibrium systems, our results are yet another 
demonstration of the important role that breakdown of spatial parity invariance may play 
in determining the universal properties. Effects of parity breakdown in pure nonequilibrium 



systems have already been elucidated in several examples in Refs. [ll|, ll2|, where breakdown 
of microscopic parity invariance is introduced by means of a reflection invariance breaking 
stochastic noise correlator. In each of those cases, universal properties depend explicitly on 
a model parameter analogous to A^x here. Noting that the disorder variance in the present 
study ([2]) is symmetric under the exchange of our model may be generalised further 

by allowing for terms in the disorder variance that are antisymmetric under exchange of i,j, 
which may separately have even and odd parity parts. We expect the coefficient of the new 
odd parity part should also appear as a tuning parameter in the critical scaling exponents. 
Calculations analogous to the above may be performed for detailed study. Lastly, it will be 
interesting to see how the predictions of Ref. jsS] for the disordered classical XY model in 
2d are modified for the cubic anisotropic 0(2) model at 2d with an impurity distribution as 
here. In addition, effects of quenched chiral disorder of the type discussed here on systems 
similar to Ref.jsJ] should be investigated. We hope our studies here will stimulate further 



theoretical and experimental studies. 
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